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Abstract 

We obtain the Bethe Ansatz equations for the broken Z^-symmetric model 
by constructing a functional relation of the transfer matrix of L-operators. This 
model is an elliptic off-critical extension of the Fateev-Zamolodchikov model. 
We calculate the free energy of this model on the basis of the string hypothesis. 
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1 Introduction 

In the two dimensional solvable lattice models with Ising-like edge interaction, the 
star-triangle relation 

pW(a, b\v, w) W(a, c\u, w) W(b, c\u, v) 

= ^2W(a,d\u,v)W(d,b\u,w)W(d,c\v,w) (1.1) 
d 
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p = p(u, v, w) independent of a, b and c 

plays a central role. In (|1.1|) , the summation on d is taken over all local states. These 
are the relations among the two Boltzmann weights W(a,b\u,v) and W(a,b\u,v). 
They live on the edges in two different directions of the two dimensional planar 
lattice. The local state variables a and b live on the sites. We denote the spectral 
parameters by u and v. 

{Fig. 1 and Fig. 2) 

Since Fateev and Zamolodchikov (T| obtained an iV-state generalization of the 
critical Ising model as a solution of the star-triangle relation (STR), there have been 
known two different off-critical extensions of this model. One is the chiral Potts 
model |J and the other is the broken Z^r-symmetric model. Both are Ising- type edge 
interaction models. The STR for the chiral Potts model was proved in |J [Q . Though 
this model has been still under investigation in || -||10||, the lack of a difference- 
variable parameterization in this model causes difficulties in analysis. Kashiwara and 



Miwa JTTJ proposed the broken Z^-symmetric model, and Hasegawa and Yamada |L2 
proved the STR for this model. Unfortunately the proof in JO] was wrong because 
of the incorrectness of the "ICU lemma" in their paper. 

In this Paper, we study the eigenvalues of the transfer matrix v, w) of the 
broken Z^-symmetric model, 

M-l 

Hu,v,w)^ a M M _\ = II W(b j ,a j \v-w)W(a j ,b j+1 \u-w), (1.2) 

i=o 

and calculate the free energy of this model. 

(Fig. 3) 

The local state variables take their values in Z/NZ. Throughout the Paper, we deal 
with the case of iV odd, N = 2n + 1. The Z2-symmetry of the Boltzmann weights 

W(a,b\u) = W(N -a,N -b\u), W(a, b\u) = W(N -a,N - b\u), (1.3) 



2 



ensures that the eigenvalue r = ±1 of the spin reversal operator 1Z is a good quantum 
number, where 1Z G End ((C )® M ) is defined by 

M times 

K = R®R®---®R, R vf ] = vfflj , (1.4) 

which satisfies 7Z 2 = 1. The vectors Vj (j G Z/iVZ) constitute an orthonormal basis 
in C^. In the homogeneous case u = v, we show first that any eigenvalue f(u) of 
$(u) = $(u, u, 0) can be written as 

M 2nM L 

j>{u)p{\ — u) 

n n 1 

p(u) = Y[e 2 (u-(2j-l) v \T/2), v = -, A = - - 77. (1-6) 

See Appendix A for the notation of the theta functions. The zeros {u%, ■ • • , u 2n M} of 
y?(-u) are described as follows 

2M 2nM, 



„(») ^ f ;' 0) p< A > , ) n Trr j ';r . n- 5 * 



#1 (u — Uj 


r/2) 




r/2) 



gifa + A/2]r/2) 
Mv h -\/2\T/2) i 



di{v k 


- Vj + r)\T/2) 


8i{v k 


- Vj -r)\T/2) 



(1.7) 



u fc = u*-^ for fc = 1, • • • , 2nM, (1.8) 
= — mod(Z©-Z). (1.9) 



We call the equation (|1.7| ) the Bethe Ansatz equation. The condition (|1.9|) follows 
from the double periodicity of ip(u) discussed in Section 4. We obtain the Bethe 
Ansatz equations above through a functional relation ( 1.1 3| ) for the transfer matrix 



of L-operators. These L-operators L(u) G End(C Ar (g> C 2 ) were originally constructed 



by Sklyanin [Tj| |1J] as a solution to the relation, 



L 01 (u - v) L 02 (u - w) RlUv - w) 
= Rl 2 / (v-w)L 02 (u-w)L 01 (u-v) on ® C 2 <g> C 2 , (1.10) 

(Fz<?. 4) 

where the upper indices 0, 1 and 2 mean that L^(u) acts only on the z-th and j-th 
components of ® C 2 ® C 2 and as identity on the other component. We denote 



the i?-matrix of the eight- vertex model by Rsv(u) P^|[I^|. We consider the transfer 
matrix C(u) of these L-operators, 



C(u) = tr C 2(L 0M (u)L 1M (u) ■ ■■L M ~ 1M (u)), (1.11) 

M times 



L 0M {u)L 1M {u)---L M ~ lM {u) e End(C N ® ■■■(g) C^C 2 ). (1.12) 

We derive the functional relation 

C(X — u — 1/4) = C(u) M (^f(u) M <S>(u -T]) + (-/(A - m)) M $(m + rj^j , (1.13) 

/M = ft(2 M | T ) 9 ''"V'lf , (1.14) 
rv ^ m fl fl i/ n . ~ ~ iMl/g) frfo + (2j - lMr/2) 

c(«) = [^4](oir)n , l(u+2i , |T/2) , 2(M _ (2j _ 1) ,| T/2) . d- 15 ) 

by the method which Baxter employed to solve the eight- vertex model |T5| jn| [fP7| . 
The functional relations corresponding in the chiral Potts model and in the RSOS 
model associated with the eight- vertex were obtained in |T3) and [|T^], respectively. 



We calculate the free energy of the broken Z^-symmetric model under the hy- 
pothesis that the solution of the Bethe Ansatz equations corresponding to the ground 
state consists of "strings of length N — 1". We show that, in the infinite lattice limit, 
the centers of these strings are distributed on the imaginary axis with the density 

9 s (2y/-LNw\NT 



« m) = 2N ^ mT K(2j=iN„\N T y -3 £w< 3- (116) 

where r = \J — Ik, and the free energy per site is 

oo sinhf u] sinhf ( u)) sinhf— ^-n) 

F{u) = _ E ^^J — u w ; u ^i . en?) 



'=i /cosh(-) cosh 2 (-^) 



The last expression agrees with the result of Jimbo, Miwa and Okado obtained 
by the use of the inversion-trick, and in the trigonometric limit of k — > oo it recovers 



the result of Fateev and Zamolodchikov JTJ and Albertini [10]. 

The organization of this Paper is as follows. In Section 2, we review necessary 
facts about the i?-matrix of the eight-vertex model, Sklyanin's L-operator and the 



broken Z^-symmetric model. We derive the functional relation ( |1.13| ) in Section 3. 



After showing commutation relations among $(«), C(v ) and 1Z, we obtain the Bethe 
Ansatz equations in Section 4. We calculate the free energy of the model under 
the string hypothesis in Section 5. Finally, in Section 6, we conclude with a brief 
discussion. We fix the notation and list the formulae for theta functions in Appendix 
[A[ Miscellaneous properties of the Boltzmann weights are summarized in Appendix 
B. We devote Appendix C to the proof of the commutativity between $ and C. 

2 Review of the Broken Zjy Symmetric Model 

j 

We fix the notation for matrices. We denote the vector *(0, • ■ • , 1, ■ • • , 0) in C m by 
v ( m \ j = o, 1,.. - ,771-1, and the matrix elements of A G End(C mi <g>C m2 <g> ■ • ■ <g>C m; ) 
by 

mi-lm2-l J7ij— -1 
31=0 j 2 =0 ji=0 



In [[L3[|[I1||, Sklyanin constructed the L-operators L(u) G End(C <8> C ) for the eight 



vertex model satisfying fll.10 ). The R- matrix of the eight- vertex model Rsy(u) |15| ]T6| 



is a solution of the Yang-Baxter equation, 

R%(u - v) R° 8 2 v (u - w) RlUv - w) 
= R\l{v-w)R%{u-w)R%{u-v) on C 2 ® C 2 ® C 2 . (2.1) 

Its non-zero matrix elements are 

r° \(u) = r\° q (u) = [e^^e^ie^iu + v), 
R%(u) = r}°(u) = [e^iv^e^u^e^u + v), 

R™(u) = RlKu) = [drfMvWtfMuWtfMu + T]). 
{Fig. 5) 

The other elements not specified above are all zero. Here we denote 9i{u)6±(u) by 

Ti 

[6 , i6 , 4 ](m) for short. We usually suppress the elliptic modulus r. When 7] = — , 



N = 2n + 1, the L-operators L(u) in ( |1.10| ) have a "cyclic" representation. In this 
representation, L(u) factorizes elementwise as 

L%{u) = K l b a {u)K^ a {u) ) (2.2) 

where 



K\ a +(T {u) = (-iy^ 2 [e 1+j e^j}(u - aarj + J), (2.3) 
^» = Gr 1 ^ 1 ^!*), G a = (^y^) 1/2 , (2-4) 

for j — 0, 1, a, b — 0, 1, • • • , N — 1 and cr = ±1. 



6) 

The factors K^{u) and K^ b a (u) are zero unless |a — fo| = 1. We can identify these 



i^(u)'s as the intertwining vectors appearing in the vertex-face correspondence [21| 
-|p6|| . Even in the Fateev-Zamolodchikov model these K(u) y s are different from the 
3-spin object V's in [0]@ by definition. Their Vs are defined by the Fourier trans- 
formed images of the product of two Boltzmann weights. These two objects, however, 
should have an intimate relationship, because the transfer matrix £ in the Fateev- 
Zamolodchikov model is also constructed from V's 0. 

Under the Z 2 -transformation which sends a to N — a, they change as 

K jN N Z a b {u) = (-iy^ +1 K ja b (u), K»Zt{u) = (-ly^^K^u). (2.5) 



They satisfy the unitarity relations 
i 



Y,G a K 3 b a {u + X)K^ c {u) = 6 ac [9 2 9 3 e 4 \(0)G b 9 2 (2u), (2.6) 

3=0 
N-l 

Y^GaK^u + ^K^u) = 8 i5 [e 2 9^mG b e 2 {2u), (2.7) 

a=Q 
N-l 

Y,G b K>(u)K»( u + \) = 6 ij [e i e 3 e 4 \{o)G a e 2 (2u). (2.8) 

6=0 

{Fig. 7) 
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In ||12|| , we determined the Boltzmann weights W and W of the broken Zjv- symmetric 
model by the relations 



W(a, b\u, v) ]T K*{u - w) Kf b (v - w) 



j=0 



^K^ a (v-w)K j d b (v-w)W(c,d\ 



u,v) 



(2.9) 



j=0 



N-l 



J2 W(a, b\u, v) K^u - w) K ]C b \v - w) 



b=0 



N-l 



K^v - w) K> b a {u - w) W(b, c\u, v) 

6=0 

(Fig. 8 and Fig.9) 



(2-10) 



From the above relations and fl2.5| ), we have the Z 2 -symmetry 
W(a, b\u, v) = W(N -a,N- b\u, v), W(a, b\u, v) = W(N -a,N- b\u, v). (2.11) 



The equation ( |2.9| ) implies that W(a, b\u, v )=W(a, b\u — v) and that 

W(a + l,b+l\u) 6 3 (u + (a + b + l)r]) 

W{a,b\u) 6 3 (u- (a + b+l)rjy 

W(a + l,b-l\u) 6 2 {u + (a - b + 1)77 



(2.12) 



W(a,b\u) 



9 2 (u - (a -6 + 1)77)' 



The crossing symmetry 



W(a, b\u) = G a G b W(a, b\X - u) 



(2.13) 



holds in this model. In the paper |S7|], Hasegawa proved the crossing symmetry only 
from ( |2.9| ) and the unitarity relations ( |2.7| ) and ( |2.8| ). We thus have 

W{a + l,b + l\u) _ G a+1 G b+1 6i(u - (a + b)rj) _ 

(2.14) 



W(a,b\u) G a G b 6 4 (u + (a + b + 2)?])' 

W(a + l,b-l\u) G a+1 G b ^ 9 1 (u - (a - b)rj) 



W(a,b\u) G a G b 6 1 (u + (a-b + 2)r ] y 

without directly solving ( |2.10|) . We can see from ( p.!2| ) and ( |2.14| ) that the Boltzmann 
weights satisfy the reflection symmetry 



W(a, b\u) = W(b, a\u), W(a, b\u) = W(b, a\u). 
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Defining T^\a\u) and T ( k '(a\u) by 



(-)/ 



T ( k + \a\u) 



n 



and T 



(-). 



am 



Tl +) (a|A- M ) 



(2.15) 



=i # fe ( u _ ( 2 j _ 1)^) 

respectively, the solutions to the recursion relations ( |2.12| ) and ( |2. 14j ) under the nor- 
malization W(0, 0\u) = W(0,0\u) = 1 are 



W(2a, 2b\u) 
W{2a,2b\u) 



Tt\a-b\u)Tt\a + b\u), 

G 2a G 2b Tt\a - b\u) Tt\a + b\u). 



Here all local state variables are to be read modulo N. See Appendix || for details. 



Hasegawa and Yamada in 0] established the star-triangle relation (STR) in the 
broken Z Ar-symmetric model, 



p W(a, b\v — w) W(a, c\u — w) W(b, c\u — v) 



N-l 



]T W{a, d\u - v) W(d, b\u-w) W{d, c\v - w), 

d=0 



(2.16) 



where p is a scalar function independent of a, b and c. 



3 Functional Relation 

In this section, we consider the transfer matrix of the L-operators and construct a 
functional relation for it. In the course of the calculation, we utilize the factorization 
property of L into K J s ([2.2|). We define a 2-by-2 matrix L(a, b, \u, v, w) by 

/ 



L(a, b, \u, v, w) 



K b a (u - w) K 0b a (v - w) K b a (u - w) K u a (v - w) 
\ K b a {u - w) K 0b a (v - w) K b a {u - w) K lb (v - w) 

Then the transfer matrix C(u, v, w) of L-operators on the lattice of width M with the 
periodic boundary condition is 

c{u,v,wt b -:^ 

= tr^L(a ,bo\u, v,w)L(ai,b 1 \u,v,w) ■ ■ ■ L(clm-i, &m-iK v, w) j (3.1) 

M-l 

= E II K ij+ X{u-w)K i %{v-w). 



(Fig. 10) 

The final goal of this section is to establish the functional relation 



£(\ - u, A - v, w + 1/4) $(w, v, w) 



( 



C(u, v, w 



f(u, V, w) $(«, V, W + 7]) 



y +(—f(X — V, \ — U, —w)) &(u,V,W — T]) 

where we define C(u,v,w) and f(u,v,w) by 



, (3.2) 



C(u,v,w) = [6 2 6 3 d4(0) [Tt ] Ti +) ]{n\u - w) [T^T^Mv - w), (3.3) 

[9i9 4 ](v- w + i]) 



f(u,v,w) = 0i(2u-2w) 



(3.4) 



[e 2 e 3 ](u-w) 

This functional relation reduces to ( p.. 13 ), ( 1.14 ) and (|1.15|) in the homogeneous case, 
u = v. We achieve this goal by the method a la Baxter ||i~5| flBfl [ft"7ll , which states the 
following: Suppose that we can find C-valued functions (fi^\b\u,v,w) {y = 0,1,2,3 
and b G Z/NZ) and matrices Pj G End (C 2 ) for j G Z/MZ, which satisfy 

AT-l 
6=0 

(jjp (a\u,v,w) <pf\a\u,v,w) 

i <f)j\a\u, v,w) 

Then we have 



for a G Z/NZ and j G Z/MZ. (3.5) 



bo,---,bM-l 



M-l M-l 

] { (j)j(aj\u,v,w) + Y[ <pf\aj\u,v,w). 



(3.6) 



i=o 3=0 
Defining vectors ifjM(u,v,w) G (C N )® M by 

M-l 



ip( u \u, v, w) 



J=0 



(3.7) 



we can write (13.61) as 



£(w, f , w) if)(°\u, v, w) = if}^\u, v, w) + ip^iu, v, w). 



(3.8) 
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In the following, we will find a family of solutions to ( |3.5| ) 

<$\b\u,v,w) = 0W(c^,6,c i+ i|u,i;,«;) (u = 0,1,2,3 and b G Z/NZ), 

P, = P( Ci ), 

labeled by j (c , Ci, • • ■ , cm-i) \ Cj G Z/NZ for j G Z/MZ}. This gives rise to 
N M vectors ip("' labeled as ip^(u, v, vj) c °' c1 ''"' Cm - 1 . We will also prove that ip's are 
proportional to the row vectors of the diagonal-to-diagonal transfer matrix $ of the 
broken Z^-symmetric model, 

^°)(A — u,X — v,w + mZZ'Zzl = $ K w)ZCZZ\, (3-9) 
^ (1) (A A -«,«; + 1/4) 

= (C(u, v, w)f(u, v, w)) M <S>(u, v,w + v)tZ'Z~-l> ( 3 - 10 ) 

^(X-u,X-X-v,w + l/4)^^Z\ 
= (-C(u, v, w)f(X -v,X — u, -w)) M <l>(u, v,w- ri)Z;:.Zz\. (3.11) 



The results ( |3.8| ) to (|3.11|) altogether implies the functional relation (|3.2| 
Now we start to solve (|3~5|) . We write the matrix elements of Pj as 



Pj 



(0) (2) 
Pj Pj 

(1) (3) 
Pj Pj 



and its first column vector t (pf\pf^) as pj. Multiplying Pj to ( |3.5| ) from the left and 
taking its first column, we have 

N-l 

<jyf'{b\u, v, w)L(a, b\u, v, w) p 7 -+i = (pp(a\u, v, w) p, for a G Z/NZ. (3.12) 

6=0 

For later use, we define the functions A*(±), A? ± \, 5* and 5* by 

A* (±) (p, a \u)= pW i^/t 1 («) " ^o a f («), 
A^ ±) (a,p|w)= -^ 0a a 1 («) + K la ^ 1 (u) , 
*.(a|u) = K Q a -\u) Kffiu) - ^» K 1 a -\u), 
5*(a\u) = K Qa - a l (u) K^+^u) - K QaJ ^{u) K la ~ l (u). 
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The equations ( |3.12| ) constitute a system of 2N homogeneous linear equations in 
4>f\a\u,v,w) and (j)j(a\u,v,w) for a G Z/iVZ. It has a non-trivial solution if and 
only if the determinant of its coefficient matrix vanishes. Demanding this condition, 
we obtain 

N-1 

Yl A*(_)(p i; a|w - w) A*_ } (a,p j+1 \v - w) 

a=0 

N-1 

+ JjA„,( + )(p J -,o|u-«;)A* + j(o,Pj + i|i;-«;) = 0. (3.13) 

a=0 

Later we find that the equation ( |3.13|) restricts p^- to a discrete set of values. When 
we parameterize Pj as 

. I p<°>(c)\ ( [Ma] (^) ) 

P 3 = P(cj), P(c) = = , (3.14) 

and denote the dependence on p(c) simply by c and U a = u + ar], A's and S's become 

A» (±) (c,o|u-l/4) = [ ^ ](0) 9 2 (U ±CTa )6 3 (U TCTa ), 
A* {±) (a,c\u-1/A) = [8 2 e 3 ](0)9 1 (U ±aTc+1 )e 4 (U ±a±c+1 ), 

S*(a\u-l/4) = [W 4 ](0)G^ 1 (2[/ 1 ). 



Under the parameterization ([3.14 ), the condition ( |3.13|) holds if and only if either c^-'s 



are all integers, or all half-integers. We restrict ourselves to the case that Cj s are 
all integers, because only in this case the relation Q3.8Q gives the functional relation 
fl3.2|). The system of equations ( |3.12| ) involves not all 0's but only <f>f\a + l\u,v,w), 
(a— l\u, v, w) and <Pj(a\u,v,w). Expressing <f>j(a + l\u,v,w) and <pp{a\u, v,w) 
in terms of 4>f (a — l\u,v,w), we obtain 

4>f\a + l\u,v,w) A* ( _)(c i ,a|w - w) A*_ } (a, c j+1 \v - w) 



<f)f\a - l\u, v, w) A *(+)(cj, a\u-w) A* (+) (a, c j+1 \v - w) ' 



(3.15) 



^ 1)(a|w) -6Mu- W ) A I- ){ " ,Cj+1 ! V ~ W l (3-16) 



<j>f\a-l\u,v,w) '" ' A*(+)(cj,a\u- w) 

From (|3.15| ) and (|3.16|) , we can write (f)j(a\u,v,w) [y = 0, 1) as 

<j)j(a\u, v, w) = (j)( u \cj, a, Cj + \\u, v, w), (3-17) 
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where the function (j)^ v '(c, a, c'\u, v, w) is independent of j. Taking the determinant of 
both sides of ( |3.5|) , we find 

det(P j+1 ) 



(5*(cj, a\u — w) 5* (a, Cj + i\v — w) 

ut '' l 1 j i 

11 11 11A 

(3.18) 



det(P,) 

4>^(cj, a, Cj + \\u, v, w) (f>j \a\u, v, w) 



1, c j+1 \u, v, w) (o) (cj, a + 1, Cj + x\u, v, w) 



Cj, a 



We set det(Pj) to unity without loss of generality. Then the equations ( |3.15| ), ( |3.16| ) 
and ( |3.18| ) give 

^\ Cj ,a,c j+1 \u,v,w) __ 6 , {alv _ w) Au(c j ,a\u-w) (3 1Q) 



0(°)(cj, a — 1, Cj + i\u, v, w) A*? + ,(a, Cj + %\v — w) ' 

where we write <pf (a\u, v, w) as (j0{cj, a, Cj + i\u, v , w). The relations Q3.15|) , ( |3.16| ) 



and (|3.19|) recursively determine (fy£> 's. We abbreviate u — w + ar/ and v — w + 777 to 
A a and _B 7 respectively. Comparing ( |3.15|) with (|2.12j) and (|2.14 ), we have 



(f>W (a, b + 1, c\u, v, w + 1/4) W{a, b + l\A ) W(b + 1, c\B ) 



(3.20) 



(o) (a, b, c\u, v,w + 1/4) = W{a, b\A Q ) W(b, c\B ). (3.21) 



(a, 6-1, c\u, v,w + 1/4) W{a, b - 1\A ) W(b - 1, c\B ) ' 
Hence we find that (p^ is a product of the two Boltzmann weights, 



From (|3.15| ) and (|3.16|) , we obtain 



(1) (a,6 + l,c\u, v,w + 1/4) _ (o) (a,6 + l,c\u - r],v - r],w + 1/4) 
0W(a, 6 — 1, c\u, v , w + 1/4) 0(°)(a, 6 — 1, c\u — 77, v — 77, w + 1/4) 

The same procedure for <// 2 ) yields 

(2) (a, 6+1, c|w, u, w + 1/4) _ (o) (a, 6 + 1, c|w + 77, v + 7/, w + 1/4) 
0( 2 )(a, 6 — 1, c|w, f , w + 1/4) 0(°)(a, 6 — 1, c|w + 77, t> + 77, w + 1/4) 

By ( ET21D , (|02D and ( pg ), we can write and (2) as 



(3.22) 



(3.23) 



</>M{a,b,c\u,v,w + l/4) = f ac {u,v,w)W{a,b\A_ 1 )W(b,c\B. 1 ), (3.24) 
(p (2) (a,b,c\u,v,w + 1/4) = ^(u, u, «;) W(a, 6|Ai) F(6, c|5i), (3.25) 



where / ac and g ac are functions independent of 6. The equations ( |3.16| ), ( |3.21| ) and 
( |3.24j ) determine f ac (u,v,w) as 

f ac (u,v,w) 
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[W 4 ](0) 6 1 (2A ) 
C(u,v,w)6 1 (2A ) 



G b giCg-6+c+i) g4Cg-fe-c+i) g^a, & - 1^0) - 1, 4go) 
Gt-i e 2 (A a _ b )0 3 (^-a-b) W(a, &|A_ X ) c|S_i) 

- C{u,v,w) f{u,v,w), 



[02e 3 ](A o ) 

where C(u,v,w) and f(u,v,w) were given in ( |3.3|) and (|3.4j ), respectively. The last 
equality is due to (|B.4|) . In the same way, we obtain 



g ac (u, v, w) = -C(u, v, w) f(X -v,X-u, -w). 



The equations (]3.24j) and Q3.25|) become 

(1) (a, b, c\u, v,w + 1/4) = C{u, v, w) f{u, v, w) W{a, 6|A_i) W{b, c|S_i), (3.26) 

(2) (a, 6, c|u, v,w + 1/4) = f, w) f(X — v,X — u, -w) W(a } b\Ai) W(b, c|Si). 

(3.27) 

Substituting ( ggg ), (f£2§) and ( gg7| ) into the definition Q of «, w) and 

using the crossing symmetry ( |2.13|) , we obtain (|3.9p, ( 3.1UQ and ( |3.11| ). We have the 
functional relation (13. 21) result. 



4 Bethe Ansatz Equations 

In this section, we give commutation relations among 7Z, C(u) and <&(v), and reduce 
the functional relation (|3.2|) to the functional equation among their eigenvalues. After 
discussing some properties about the zeros and poles of the eigenvalues of $(«), we 
derive the Bethe Ansatz equations ( |1.7| ) , (Of) and ( |1.9| ) for the broken Z Ar-symmetric 
model. 

First we have 

C(u, v, w') $(it, u, w) = $(it, v, w) C(v, u, w'). (4.1) 

We give a proof in Appendix O. In the case of the homogeneous system, i.e., u = v 
and w = w', the equation ( |4.1| ) means the commutativity of two transfer matrices 



C(u) = C(u,u,w) and $(it) = u, w), 



C(u) , J = 0. (4.2) 
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The star-triangle relation ( |2.16| ) gives 



, = 0, 



(4.3) 



and the LLR = RLL relation ( |1.10| ) guarantees 



Ciu) , C(v) = 



(4.4) 



The relations (fO|), ( |4.3j ) and ( ^.4| ) make it possible to diagonalize £(w) and $(t>) si- 
multaneously by eigenvectors independent of the spectral parameters u and v. Fixing 
one of the eigenvectors and denoting the corresponding eigenvalues of C(u) and $(t> ) 
by l(u) and tp(v) respectively, we can rewrite the functional relation (|3.2|) as 



J(A - u - 1/4) tp(u) = C(u) M ^f(u) M <p(u -rj) + (-l) M f(X - u) M V (u + rj) j , (4.5) 
where from ( |3.4| ) and ( |3.3| ), f{u) and C(m) are 

f(u) = O^u) ^l^^ and C{u) = [6 2 d z dm[Tt ) Tt ) Tt ) Tt\n\u). 



(4.6) 



The next step to derive the Bethe Ansatz equations is to examine the quasi-periodicity 
property of y?(w)- We have the following relations 



$(tt + 1) = 



7\ 



(4.7) 
(4.8) 



We have defined 7?. in (1A). The periodicity (17) is obvious from the definition (1.2) 
of <£> and the periodicity ( P-5|) of W and W 7 . The other periodicity ( f4.8|) follows from 
( B.6|) . We also have the commutativity 



K , C{u) 



0. 



(4.9) 



from the Z2-symmetry of K's ( |2.5|) and the definitions (|1.4|) and ( [3.2|) . Diagonalizing 
7?. and $(tt) simultaneously with C(v), the equations ( |4.7|) , (|4.8|) and (|1.4|) give 



<p(u + l) = <p(u), <p(u + ~)=r<p(u), r = ±l, 



(4.10) 
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where r is an eigenvalue of 1Z. The poles of ip(u) are coming from only those of the 
matrix elements of $(«). We define 

P(«) = (S^YfllW(u-(23-l)v\T) 



= Y[e 2 (u-(2j-l) V \r/2), (4.11) 

i=i 

which contains all possible poles of W(a, b\u). The same does p(X — u) for W by the 
crossing symmetry Hence the set of zeros of [p{u)p{X — contains all poles of 

(p(u). By the Lemma in Appendix [A| and the double periodicity ( |4.10| ) of <f(u), we 
can write <^(u) as 

2nM 

nfii(«-«i|r/2) 

y>(u) = (consi) i=- 

p(A — u) 



2nM I - r T 

Uj = nMX + — — mod L © —Zj). 

i=i ' 4 2 



The initial condition $(0) = Id ( |B.3j ) determines the normalization of const, 

v{u) = ( mm r^tyy), ( 4.i2) 



j?(tt)p(A — w) 



5=1 



Assuming C(uj) ^ in (|4.6|) and substituting u fc (fc = 1, ---,2nM) into (|4.5|) , we 
have 

f(u k ) M (p(u k -r,) + (-1) M /(A - « t ) A VK + »7) = 

for fc = l,-..,2nM. (4.13) 

We further assume that Uk (k = 1, • • • ,2nM) are neither zeros nor poles of f(u), 
f(X — u) and tp(u ± 77). Then ( [4.13 ) becomes 



f(X-Uk)p(u k -r])p(X-Uk + ri) J ' ~ =1 0i(ufc - uj - t]\t /2) 



for fe = 1,--.,2tiM. (4.14) 



We can write /(tt) in ( P~5D by ( [OP as 



f { "^ = Z^ e Mr/2) e l (u + V \r/2). 
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Then the left-hand side of (|4.14|) reduces to 
by A/2, i.e., putting 



' giMr/2) 
9 1 (u k -X\r/2) t 



2AI 



. After shifting u k 



Vk = Uk — — for k — 1, • • • , 2nM, 



we obtain the Bethe Ansatz equations for the broken Z^-symmetric model, 



Mv k -X/2\r/2) 



2M 



2nM 

E 



V; 



jA0i(vk-Vj-v\r/2) 

for k = l,---,2nM, (4.15) 

mod(Z©^Z). (4.16) 



5 Density Function and Free Energy 



In this section, we will calculate the free energy of the broken Z^r-symmetric model 
from the Bethe Ansatz equations under the three assumptions concerning the ground 
state. One is the String Hypothesis below, and the others are about the distribution 
of string centers v a = \/—lw a and the corresponding quantum numbers I a . We 
restrict the spectral parameter u to the region [0, 1/2N\, in which all the Boltzmann 
weights are real and positive, and r to a pure imaginary number, r = ^/—Ik with k 
real and positive. 

By a string of length I and parity v (= or 1) with its center v a , we mean the 
following set, 



v ad = v a + (2j -1-1)^ + ^ mod(Z©^Z) 



(5.1) 



for j — 1, 2, ■ ■ • , I, and v a : pure imaginary 

We suppose that the following hypothesis holds in the infinite lattice limit |10[ [p8 
String Hypothesis for the ground state 
The solution of the BAE's ( \4-lo] ), {vj,j = 1, • ■ ■ , 2nM}, cor- 
responding to the ground state consists of strings of length 
1 _ 

N — 1 and parity . 
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More precisely, for finite systems the solutions of the BAE's may have deviations 
from strings. The hypothesis asserts that these deviations vanish in the infinite lattice 
limit. In the course of the following calculation we deal with the solutions of the 
BAE's as if they were genuine strings, since we are interested in thermodynamic 
quantities. Because all the matrix elements of $ are real and positive, the Perron- 
Frobenius theorem 2^1 shows that the ground state belongs to the sector of zero 



quasi-momentum []10|1 - The hypothesis implies that the ground state also belongs to 
the sector r = 1, and that the corresponding solutions are made up of M strings of 
length 2n. We denote them by 

i - f-rp +1 



v aJ = V^lw a + (iV-2j)| + 

for a 



mod (Z © -Z) 



M and j — 1, • • • , 2n, 



where wJs are all real and taken as 



-7 < Wi < w 2 < ■ ■ ■ < % < -• 



Then the BAE's for the ground state becomes 

2M 



6^v p>k + \/2\T/2) 
Mvm-W\t/2). 



' iw+infT gl( ^' fc ~^' i + r?|r/2) 

for /3 = 1,-..,M, k = l,---,2n. 



(5.2) 



Multiplying (|5.2|) over k = 1, • • • , 2n, we have 



2n 



\k=l 



n 

2N { 



l-(-l) 



n+l 



2M 



2n-l 



a=l 



where a) is 



2n y 




4 






n 






- w«, 




)U 


[wp - 


0i (a 




—lw 


r/2) 


01 (a 


+ V 7 


—lw 


r/2) 



2r^ 



Taking the logarithm of (|5.3| ) and dividing it by 2a/— IttM, we have 



(5.3) 



T(w p ) = ^ for /3 = 1, 



,M, 



(5.4) 
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where the quantum number Ip's are integers and 

I M 

T{w) = Ti(w) - — ^T 2 (w- w a ), 



a=l 



i , i - f-ir +l 



2n / n 1 
T x (w) = Yt(w, — (n-2k + —, 

K J ^ V ' 2iV v 2n' 4 

2n n 2n 2 
T 2 H = *( W; 0) + 2$>(™,-7)-*K— ), 

3=1 

t(tu,a) = * logx(w,a). 

-J.7T 



Albertini et al. numerically investigated the 3-state Fateev-Zamolodchikov model 0. 
Their results indicate that the String Hypothesis holds. We assume that the centers 
w of the strings are distributed densely on the interval [ — k/4, k/4] in the limit of M 
large, and that the quantum numbers Ip satisfy 

Ip +1 =I p + \ for /3 = 1, 2, • • • , 2nM. (5.5) 

These are the second and the third assumptions we make. The results by Albertini 
et al. also supports them. We furthermore conjecture that 

w a = -w M - a+1 (5.6) 

holds exactly for the ground state even in the finite lattice. This conjecture is con- 
sistent with their results. If Q5.6|) is true, we can show that 



T(«/4) - T(-k/4) = 1, (5.7) 

and this implies that 2M integers Ip's must fill the interval [ — M, M) without jumps 
if all Ip's are different. This also supports our assumption. But we do not use the 
conjecture ( |5.6| ) in this Paper. 

We now proceed to the calculation. We define the density function for w's by 

^ ^ M-*oo M(wp + i — Wp) ' ^ 



which is positive and for any integrable function f{x 

M 



J™,, 77 fM = l*J( w )p( w ) dw ( 5 - 9 ) 

a=l 4 



holds. Considering the difference of (|5.4|) for Ip + i and Ip 



M(wp +1 - wp) M(wp +1 - w p ) 
and letting M — > oo, we obtain 



Wp + i - Wp 



, . dT(w) dTAw) 1 [j dT 2 (w — w) ,_. ,_ 
p(w) = — = — / p(w) dw. 



dw dw 2 J 

By (|A.l ) and (|A.4 ), we can expand t(w,a) as 



dw 



(5.10) 



t(w, a) 



8n 



4ir 



{a}w 



(2{a} 



1 , h 
log — 



{a} H — w 

K K 



-1 



0i 



-{a} w 



oo sinf w) sinhf (\a\ 



k sinhf ) 



We denote the fractional part of x by {x} = x — [x], [x] being the Gauss symbol. 

dTAw) dT 2 (w) 

When we write the Fourier expansions of — , — and p(w) as 



dTj(w) 



fc=— oo 
oo 



dw 
p(w) 

the integral equation ( |5.10| ) gives 



A jk exp 



dw dw 



K 



-w^j for j = 1, 2, 



^ /4V-17TA; 
2^ Pk exp(^ u> 

fe=— oo 



K 



Pk 



K A 

i + iM,k 



The coefficients A-^ are 



4n 



4 sinh (£^ ( AT_ 1) ) cosh (^ (Ar+i: 



K 



/27ck\ /27ck\ 
sinh coshl I 

v K ' v K ' 



1 + ~^%k 



2n 
~N 



2sinh(— (iV-1)) cosh(— 



sinhf ) 

v K ' 



k = 0, 

Mo, 

fc = 0, 
k ^ 0. 
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We obtain the density function for strings 



2 oo exp( w 



fe=— oo 



coshf^-) 



Using (|A.lj) and (|A.5|) , we can rewrite it as 



.1^1(01-^)^ 



2iV [0 2 3 ](o|iVr 



2w 


1 


K 


~ Nt 


2w 


1 


K 


~ iVr 



e 3 (2-/=TiV«;|iVr 



9 2 {2y/^lNw \Nr 
The free energy per site F(u) of the model is defined by 



(5.11) 



(5.12) 



where <p(tt) is the eigenvalue of the transfer matrix $(«) corresponding to the ground 
state. We can write (p(u) as 

( p(0)p(\) 



M 2n M 



fl, = { 7 ,}, 7 ,- = —(N-2j + — ) + - ^ , 

f3 X ' 3S ' 13 2N K J 2n 4 

by (|4.12|) , (|5.6|) and the String Hypothesis. Then the free energy per site is 

p(0)p(\) x 2n 1 M 



F{u) 



-log | 
log | 



\ In ^ M 

- — lim — log-D,- 

' / j=l a=l 



U,W ay 



s p(u)p(X — u) 



Since p(w) is an even function, it is enough to integrate the even part of log Dj(u, w). 
We hence have 

p(0)p(A) \ 



F(u) 



log | 



p{u)p{\ — u) 



1 2n i 



2n „k 



d log -D^(tt, w 
dw 



■ p( 1 \w)dw, 
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where D^(u, w) and p( ^(w) are 



3 V""' 



/w 
p(w) dw, 
-f 

Df\u,w) = (D j (u,w)D j (u,-w)^ 1/2 



A little cumbersome calculation yields 



1 2n 



3=1 

2^y_. ^ ^ = -e 3 (m) + e 4 (m), 



where 



i=i 4 



. . 27m / 1 
^(„) = _„(_-2« 

E 2 (u) - if Sinh ^") Sinh ( V ( 2A7 ~ U) ) C ° Sh O Sinh (^") 
/=1 I sinhf ) sinhf 

2ttZ n . , /2tt/ , 1 v \ /2tt/ , _ N . . /2tt/ 



oo sinhf sinhf (— — — w)) coshf— — (n + 1)) sinhf— —n) 

E 3 (u) = -4XX-1)' V — 7 - ^ ' K -ff- - KNk ' 

1=1 I sinhf ) sinhf ) 

sinhf u] sinhf (u —j) coshf— — (n + 1)) sinhf-— n) 



E i {u) = AY J 



/2ttI\ /2tyI\ / TV I 
1=1 I sinh sinh cosh 

Now F(u) has the final expression 
F(u) = E 2 (u) + E 4 (u) 

h Zcoshf^cosh 2 ^ 



It agrees with the result of Jimbo, Miwa and Okado [20| obtained by the use of the 
inversion-trick. In the trigonometric limit k — > +oo, the free energy formula ( |5.13 ) 
reduces to the integral 

sinh f Nttxu) sinhfA^7rx( uj) sinhfn7rx) 

lim F(u) = - dx %? 2AL 



K^OO 







x cosh( — — cosh 2 (^) 
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which also agrees with the results of Fateev and Zamolodchikov themselves ffl] and 
by Albertini [10|. The former was obtained by the inversion-trick, and the latter by 
the Bethe Ansatz method. 



6 Discussion 

The first main goal of this Paper is the functional relation ( |3.2|) . We obtain it as 
a functional relation for C(u). The diagonal-to-diagonal transfer matrix of 
the broken Z^-symmetric model appears naturally in this relation. We obtain the 
Boltzmann weights W and W of the broken Z jv-symmetric model in the algebraic 
way different from that of [ F2]| . We obtain W and W as the solutions to the relation 



( |3.12j ). In fl2]| , they are the solutions to the relations ( |2.9| ) and ( [2.10| ). 



The Bethe Ansatz equations ( |4.15| ) are the second goal of this Paper. The commu- 



tativity ( flT|) between C{u) and <3>(t>) is essential to get the Bethe Ansatz equations 
( |4.15| ) from the functional relation ( |3.2|) . It is notable that the unitarity relations 



( p.6| ) and (|2.7] ) guarantee this commutativity. This contrasts with the usual situation 
where the commutativity of the transfer matrices is derived from the STR or the 
LLR = RLL type relations. 

The Fateev-Zamolodchikov model is the trigonometric limit of the broken Zjy- 
symmetric model. It has the Z^-symmetry besides the Z2-symmetry. Hence 
the Z^r-charge q G {— n, —n + 1, • • • , n — 1, n} is a good quantum number, where 



exp(2v— Tvr— ) is an eigenvalue of the Z^-charge operator Q e End((C 

M times 



Qvf ' =exp [2^lnj-)vf ] for j = 0, 1, • • • , iV - 1. (6.1) 



Albertini [10] obtained the Bethe Ansatz equations and the formula for the eigenvalue 



(fFziu) of the diagonal-to-diagonal transfer matrix for the Fateev-Zamolodchikov 
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model. They are 



f ( ,,m\2M 2nM-2\q\ / . s 

s + A/2 \ =( _ 1)A/+1 n sK - vj + y) forfc = 1 2nM (6 . 2) 



M 2nM-2| 



«>«(«) = "r5 0) a'"' J "n ,£ ^. P00W = Jim p(«), (6.3) 



where s(tt) = sin(7ra) and is given in ( 4.11 ). The notations are slightly changed 



from |L0| to compare to our case. The main difference is the number of factors in 
the right-hand sides of the BAE's and in the expressions of the eigenvalues <p(u) and 
^fz(w). It is always 2nM in the broken Z^-symmetric model, and in the Fateev- 
Zamolodchikov model it is 2nM — 2\q\, which depends on the sector of the Z^-charge 
operator Q. This difference originates in that the Z^-symmetry holds only in the 
Fateev-Zamolodchikov model and that it breaks away from the criticality. The BAE's 
( |4.15| ) for the broken Z jv-symmetric model should coincide with those for the Fateev- 
Zamolodchikov model in the trigonometric limit k — > 00. We conjecture that the 



situation is the following. In the solution {vi, ■ ■ ■ , v 2n M} to the BAE's (|4.15| ) for the 
broken Z^-symmetric model, some of them diverge to ±\/— Too all in the same order 
in k when the trigonometric limit is taken. The half of them diverge to \J — loo, and 
the other half to — V— loo. The number of them are always even and between and 
2n. Let 2q be this number. Then q determines the sector of the Z^-charge operator 



in which this eigenvalue falls. In this situation, the BAE's ( |4.15| ) and the eigenvalue 



<p(u) in (|4.12|) surely become the BAE's ( |6.2| ) and <Pfz(u) in ( |6.3j ) respectively in the 
trigonometric limit. 

The free energy ( |5.13 ) agrees with the result of PD| by the inversion-trick. The 



string hypothesis for the ground state in Section 5 is consistent with their result. In 
our formulation, it is manifest that the free energy F(u) is doubly periodic in u, 

F(u + 1) = F(u + t/2) = F(u), 

from (|4.10|) , (|5.12|) and the fact that the ground state belongs to the sector of r = 1. 
This result of F(u) also gives the ground state energy of the one dimensional spin 
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chain Hamiltonian Ji corresponding to the broken Z^-symmetric model, 

log$(u) = Id + uH + 0(u 2 ). 

The Hamiltonian TC itself is modular invariant, T~C(t) = Tl( — ). We will report on 

r 

these accounts elsewhere. 



A Theta Function 



We summarize the necessary facts about theta functions in this Appendix. See PU[ 
for proofs. We define 9\(u\r) by 

oo 

0i(u|t) = 2g 1/4 sin(Tru) - 2q 2n cos(2nu) + q 4n )(l - q 2n ), 

n=l 

where q = exp(y/— Vkt). It is an odd function in u and has the quasi-periodicity 

^(u + llr) = 9 1 (u\r), 

9x(u + t\t) = — q^ 1 exp(2\/— lnu) 9\{u\r), 



and satisfies 



9i(u\t) 



- 1 ( ) exp 



T 



-Iff— W- 



(A.l) 



The other theta functions 92, 9 3 and #4 are defined by 



9 2 ( u \t) = 9x(u + l/2\r), 

9 3 (u\t) = -g 1/4 exp( v /Z l7ra)0 1 (u + 1/2 + t/2|t), 
9 a (u\t) = -v /z Tg 1/4 exp(v /z T7TM)0 1 (M + r/2|r). 

We abbreviate a product of theta functions of the same argument to, for example, 

[9 1 9 2 ]{u\t) = 9 1 {u\t) 9 2 (u\r), 
[020 3 04](O|r) = 2 (O|r) 3 (O|t)0 4 (O|t). 



In this notation, 
[9i9 4 ](u\t) 



T](t) 

V (Ar) 



9 1 {u\t/2), [9 2 9 3 }{u\t) 
0i(2u|2t), [0 3 4 ] (u\t) 



t](2t) 2 

T]{t) 

V(2r) 2 
V (4t) 



9 2 (u\t/2), 
4 (2m|2t) 



(A.2) 
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hold, where 

oo 

7Kr) = g 1/24 II (!-<?") 

n=l 

is the Dedekind eta function. The necessary addition formulae are 

[9 2 9 3 ](u\r) [0 2 3 ]Mt) ± [0i04](u|r) [0i0 4 ](«|r) = [0 2 3 ](O|r) 9 2 {u T «|r) 3 («±^ 
[0104] («|r) [0 2 3 ]Mr) ± [0 2 3 ]Hr) [0i0 4 ]Hr) = [0 2 3 ](O|r)0 1 («± ? ;|r)04(wT^ 



(A.3) 

In Section 5, we use the Fourier expansions, 

1_ ftfr + „|r) = _ 2 + 2 g sin(2^ - t/2» 

- i — v\t) k sin(7r/cr) 

3 (w|r) 1 ~ exp(v^T7rA;M) 

2^ — TZ \ — ' ( A - 5 ) 



9 4 (u\r) [9 3 6 4 ](0\r) ^ cos(irkr) 
which are valid for < lm(w) < r and v real. We are denoting the fractional part of 
by {x}. The expansion ( |A.5 ) is essentially the same as that of the Jacobi elliptic 



x 

function dn(u, k) 



The next lemma is fundamental. 

Lemma 1 Let f(u) be a meromorphic function which is not identically zero and has 
the quasi-periodicity property 

f(u + 1) =exp(-2^T7r B)f(u), 
f(u + t)= exp(-2 v ^T7r(A 1 + A 2 uj) f(u). 



Denoting the zeros and poles of f(u) by Ui, u 2 , ■ ■ ■ , u n and v±, v 2 , ■ ■ ■ , v n respectively, 
then we have 

n m i 

n-m = A 2 , -^Vj = -A 2 - A 1 + Br mod(Z©rZ), 

3=1 3=1 2 

and 

n 

Y[9i(u- Uj\r) 
f{u) = C exp( v /Z T7r(A 2 - 2B)r V_J 



H9i(u- Vj \T) 

3=1 

with C independent of u. 
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B Boltzmann Weights 

In this Appendix, we list some formulae about the Boltzmann weights. Solving the 
recursion relations (|2.12|) and ( [2.14|) under the normalization of W(0, 0\u) = W(0, 0\u) 



= 1, we have, for a, b = 0, 1, • • • , n, 

W(2a,2b\u) = W(N — 2a,N — 2b\u) = T ( 2 + \\a - b\ \u) T ( 3 + \a + b\u), 

W(2a, N - 2b\u) = W(N - 2a, 2b\u) = T ( 2 + \a + b\u) T^\\a - b\\u), 
W(2a, 2b\u) = W(N -2a, N - 2b\u) = G 2a G 2b T ( 2 ] (\a - b\ \u) T ( f\a + b\u), 
W(2a, N - 2b\u) = W(N - 2a, 2b\u) = G 2a G 2b T { f\a + b\u) T ( f\\a - b\ \u). 
Noting that T fc (fT) (a\u) in (|2.15|) satisfies 

Tt\0\u) = Tt\N\u), (B.l) 
Tt\N-a\u)=Tt\a\u), (B.2) 

we can extend the domain of the first argument of to all integers by periodicity. 
With this convention we can rewrite the above expression of the Boltzmann weights 
simply as 

W{2a,2b\u) = T 2 i+ \a-b\u)Tt\a + b\u), 
W{2a,2b\u) = G 2a G 2b T!f\a-b\u)T^\a + b\u). 

We have in particular at u=0 and A, 

W(a,b\0) =G- 1 G b - 1 W{a,b\\) = 1, G a G h W(a, b\X) = W{a, 6|0) = 5 ab . (B.3) 
When we write u + arj as U a , the next identity for 

Ti ai \n-a\u) _( 6 k (U (l+ai)/2 ) Tfc (,Tl) (n|w) 



Tl ai \a\u + a 2 rj) \0k(U^ 2a2+(1+ai)/2 ) 
gives 



W(a,b-l\u) 



W (a, b\u + arj) G b 



(01 ) 04 (^(a+6)+l) ) " [T 2 H T 3 H ] (n\u). 



(B.4) 
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7^ has the quasi-periodicity 

T^\a\u + 1) = exp(4v^T7rcTa 2 r/) T^\a\u + r) = T^\a\u) 



T fc w (a|u + r/2) = exp(-2 v /3 T7raa 2 r/) T^ ; (a|«), fc + fc' = mod 5. 
Hence we have 

W(a,b, \u+l) = W(a,b\u), W(a, b, \u + 1) = W(a, b\u), (B.5) 

and 

W(a, b\u + r/2) = exp(-4 v / ^T7r(a 2 + b 2 )r]) W(a, N - b\u), 
W(a, b\u + r/2) = exp(4 v / ^T7r(a 2 + b 2 )r]) W(a, N - b\u). 

We note 

W(a,b\u)W(b,c\u) 

= exp(-4v^T7r(a 2 - c 2 )r/) W(a, N - b\u + r/2) W(iV - 6, c|u + r/2). (B.6) 

C Commutation Relation between C and 



In this appendix, we give a proof of (|4.1| ). A graphical representation of this proof 



for the case of M = 2 is illustrated in Fig. 11. 



,l\ \ CQ--CM-1 

I ao---a M -i 



)aa—a,M-i 
bo---bM-i 



M-l 

= E ( II ^(^-1. c > - «0 ^te, - «0 ^v. ,!;;(" - «0 - «0 

By the unitarity relation ( |2.7|) , inserting 

J2G a , K^{w -w' + \)K™(w - w') 

we have 

[6 2 e 3 e 4 \(0)G Co 62(2w-2w')y 1 E Y, G « ,K i%( w - w ' + x ) 

ft O'"'>A/-l Vol 
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M-l 

x ( II W( Cj , bj\v - w) W{b^ Cj \u - w))K^,(w - w')K io b Z-_l(u - w') 

j=0 

M-l 

x ( II ^%(v - w>) An - w'))K^(v - w>). 



Successive use of (^9| ) and (|2.10|) yields 



([W 4 ](0) G C0 6 2 (2w-2w') 

E YP*> Kfi (w-w' + A) K'Z (w - w') 



ivi — z 

( II K^iv - w>) K%\(u - w'))K iM _Z Z \(v - w>) K iM -^(u - w>) 



X 

b 0"' b Af-i i'a' 
M-2 

X [ 

j=0 

M-2 

x ( II W(bj, dj\v — w) W(dj, bj + i\u — w) 

3=0 

x W(bM-i, a-M-i\v — w) W(a M -i, o!\u — w). 

By the unitarity relation ( |2.6|) , we have 

J2G ao K^(w -w' + XW'C^w - w>) 

j/_ 

[W 4 ](0) G ca 6 2 {2w-2w') 

then the above formula reduces 

M-l 



E ( II K^{v - w') K l %\{u - w') W(bj, aj \v - w) W(a Jt b j+1 \u - w 



b ---b M _i j=0 



bo---b M -i 

= (C(u,V,w')^U,V,w))Z: C a M M -\- 

Now we obtain 

v, w) C(u, v, w') = C(u, v, w') v, w). 



Acknowledgment 

The author would like to express his gratitude to Prof. M. Kashiwara, Prof. T. Miwa 
and Prof. M . Jimbo for their encouragement. Thanks are also due to Dr Y-.H. Quano 
for his careful reading of the manuscript. 



28 



References 



[1] V.A. Fateev and A.B. Zamolodchikov, Self-dual solutions of star-triangle rela- 
tions in Zjv-model, Phys. Lett. A92, pp. 37-39, 1982. 

[2] H. Au-Yang, B.M. McCoy, J.H.H. Perk, S. Tang and M. Yan, Commuting trans- 
fer matrices in the chiral Potts model: solutions of star-triangle equations with 
genus > 1, Phys. Lett. A123, pp219-223, 1987. 

[3] R.J. Baxter, J.H.H. Perk and H. Au-Yang, New solution of the star-triangle 
relations for the chiral Potts model, Phys. Lett. A128, pp. 138-142, 1988. 

[4] H. Au-Yang and J.H.H. Perk, Onsager's star-triangle equation: Master key to 
integrability, Adv. Stud, in Pure Math. 19, pp.57-94, 1989. 

[5] R.J. Baxter, V.V. Bazhanov and J.H.H. Perk, Functional relations for transfer 
matrices of the chiral Potts model, Int. J. Mod. Phys. B4, pp.803, 1989. 

[6] R.J. Baxter, Free energy of solvable chiral Potts model, J. Stat. Phys. 52, pp.639, 
1988. 

[7] G. Albertini, B.M. McCoy and J.H.H. Perk, Eigenvalue spectrum of superinte- 
grable chiral Potts model, Adv. Stud, in Pure Math. 19, pp. 1-55, 1989. 

[8] , R.J. Baxter, Chiral Potts model, eigenvalues of transfer matrix, Phys. Lett. 
A146, 110-114, 1990. 

[9] G. Albertini, S. Dasmahapatra and B.M. McCoy, Spectrum and Completeness 
of the integrable 3-state Potts model: a finite size study, Int. J. Mod. Phys. A7. 
Suppl. 1A, ppl-53, 1992. 

[10] G. Albertini, Bethe-Ansatz type equation for the Fateev-Zamolodchikov spin 
model, J. Phys: Math. Gen. 25, pp. 1799-1813, 1992. 

[11] M. Kashiwara and T. Miwa, A class of elliptic solutions to the star-triangle 
relations, Nucl. Phys. B275 [FS17], pp. 121-134, 1986. 

29 



[12] K. Hasegawa and Y. Yamada, Algebraic derivation of broken Z N symmetric 
model, Phys. Lett. A 146, pp. 387-396, 1990. 

[13] E. Sklyanin, Some algebraic structure connected with the Yang-Baxter equation, 
Funk. Anal. Appl. 16 No4, pp. 27-34, 1982. 

[14] E. Sklyanin, Some algebraic structure connected with the Yang-Baxter equation, 
representation of quantum algebra, Funk. Anal. Appl. 16 No4, pp. 34-48, 1982. 

[15] R.J. Baxter, Eight-vertex model in lattice statistics, Phys. Rev. Lett. 26, pp. 832- 
833, 1971. 

[16] R.J. Baxter, Partition function of eight-vertex model, Ann. Phys. (NY.) 70, 
pp. 193-228, 1972. 

[17] R.J. Baxter, "Exactly Solved Models in Statistical Mechanics", Academic Press, 
London, 1982. 

[18] V. V. Bazhanov and Yu. G. Stroganov, Chiral Potts model as a descendent of 
the six-vertex model, J. Stat. Phys., Vol. 59, Nos. 3/4, pp.799-817, 1990. 

[19] V. V. Bazhanov and N. Yu. Reshetikhin, Critical RSOS models and conformal 
field theory, Int. J. Mod. Phys. A, Vol 4, No. 1, pp. 115-142, 1989. 

[20] M. Jimbo, T. Miwa and M. Okado, Solvable lattice model with broken Z N sym- 
metry and Hecke's indefinite modular form, Nucl. Phys. B275 [FS17], pp. 517- 
545, 1986. 

[21] R.J. Baxter, Eight- vertex model in lattice statistics and one-dimensional 
anisotropic Heisenberg chain I, some fundamental eigenvectors, Ann. Phys. 
(NY.) 76, pp. 1-24, 1973. 

[22] R.J. Baxter, Eight-vertex model in lattice statistics and one-dimensional 
anisotropic Heisenberg chain II, equivalence to a generalized ice-type lattice 
model, Ann. Phys. (N.Y.) 76, pp.25-47, 1973. 

30 



[23] R.J. Baxter, Eight-vertex model in lattice statistics and one-dimensional 
anisotropic Heisenberg chain III, eigenvalues of transfer matrix and Hamilto- 
nian, Ann. Phys. (N.Y.) 76, pp.48-71, 1973. 

[24] V.Pasquier, Etiology of IRF model, Comm. Math. Phys. 118, pp.357-364, 1988. 

[25] K. Hasegawa, Crossing symmetry in elliptic solutions of the Yang-Baxter equa- 
tion and a new L-operator for Belavin's solution, J. Phys. A26, pp.3211, 1993. 

[26] Y.-H. Quano, "Sklyanin algebra and solvable lattice model" in the proceedings 
of the Workshop on "New Aspects of Quantum Field Theories" , S. Nojiri ed., 
INT-T-513, pp. 29-44, 1992. 

[27] K. Hasegawa, On the crossing symmetry of broken Z^-symmetric solutions of 
the Yang-Baxter equation, in "Representation theory of Lie groups and Lie al- 
gebras", the proceedings of the Kawaguchi-ko conference, T.Kawazoe, T,Oshima 
and S.Sano eds, World Scientific, Singapore, pp. 22-58, 1992. 

[28] M. Takahashi and M. Suzuki, One-dimensional anisotropic Heisenberg model at 
finite temperature, Prog. Theo. Phys. vol48, 6B, pp.2187-2209, 1972. 

[29] F.R. Gantmacher, "Matrix Theory", Chelsea Publishing Co. New York, 1959. 

[30] D. Mumford, "Tata Lectures on Theta I", Birkhauser, Boston, 1983. 

[31] E. T. Whittaker and G.N. Watson, "A Course in Modern Analysis", 4th ed., 
Cambridge Univ. Press, 1958. 



31 



Figure caption 

Figl Graphical representation of 

W(a, b\u, v) and W(a, b\u, v) 
Fig2 Graphical representation of the STR 
Fig3 Graphical representation of 3>(it, v, w) 
Fig4 Graphical representation of (1.10) 
Fig5 Graphical representation of the YBE 
Fig6 Graphical representations of 

Kfa(u — w) and K&(y — w) 
Fig7 Graphical representations of 

unitary relations between K's 
Fig8 Graphical representation of (2.9) 
Fig9 Graphical representation of (2.10) 
FiglO Graphical representations of £(u, v, w) 
Figll Graphical representation of a proof of 

C(u, v, w') <E>(w, v, w) = 3>(it, v, it;) £(u, v, w') 



32 



W(a,b\u,v) 



W(a,b\u,v) - 




Fig.l Graphical representations of 

W(a, b\u, v) and W(a, b\u, v) 
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Fig. 2 Graphical representation of the STR 
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Fig. 3 Graphical representation of <&(u,v 
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Fig. 5 Graphical Representation of the YBE 
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Fig. 6 Graphical representations of 
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= s ac [e 2 e 3 e 4 }(o)G b e 2 (2u 
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Fig. 7 Graphical representations 

of unitary relations between ICs 
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8 Graphical representation of 
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Fig. 10 Graphical representations of £(u, v, w) 
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w' 




Fig. 11 Graphical representation of a proof 
of C(u, w r ) w) = w) C(v, u, 
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